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Abstract 

The existence of the Hamiltonians of the renormalized point interactions in two and 
three dimensional Riemannian manifolds and that of a relativistic extension of this model 
in two dimensions are proven. Although it is much more difficult, the proof of existence 
of the Hamiltonian for the renormalized resolvent for the non-relativistic Lee model can 
still be given. To accomplish these results directly from the resolvent formula, we employ 
some basic tools from the semigroup theory. 

1 Introduction 

Typical field theory problems require a concept known as renormalization, which is a way 
of rendering infinite quantities to finite values to get physically meaningful results. This is 
a very hard problem, and it would be interesting to find some simple models in which the 
ideas can be tested in depth and a mathematically sound description can be presented as 
much as possible. This will illuminate the underlying mathematical and physical ideas in more 
complicated models. 

There are indeed some simple models introduced in the past. One of them is the Dirac-delta 
potentials in quantum mechanics which was first studied rigorously by Berezin and Fadeev [1] 
and later discussed extensively by Albeverio et al [21 E]. These works show that Dirac-delta 
potential can be understood from the self-adjoint extension point of view, hence all could be 
made mathematically sound. Many body version of this problem on is known as the formal 
non-relativistic limit of the scalar field theory in (2-1-1) dimensions. All these are extensively 
discussed first in the unpublished thesis of J. Hoppe [1]. A similar model is the non-relativistic 
Lee model, which exhibits an additive divergence. We are not aware of a mathematically 
rigorous discussion of this model. Physically the relativistic simplified version of the Lee model 
is more important and there is quite of a bit of work from a nonperturbative point of view to 
understand the physics behind it (see the references in [5]). The approach we follow is introduced 
in [H] by Rajeev, and recently we have introduced the generalizations of these models on to 
manifolds [6l [TJ El HI [10] . The rigorous understanding of the existence of the Hamiltonian left 
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aside in our previous works. We would like to address this issue in the present work. There is 
a general approach which is exposed in the excellent book by Albeverio and Kurasov and 
it should be applicable in the Dirac-delta functions case for the manifolds. However, we will 
use an alternative approach. The advantage we have is the following, the self-adjoint extension 
point of view becomes complicated when we discuss field theories, it is usually hard to give a 
meaning to operator valued distributions and their extension theory is even more delicate. The 
other alternative which uses resolvent convergence of cut-off Hamiltonians is problematic when 
we use other regularization schemes, e.g. the powerful dimensional regularization. This is why 
we want to utilize a direct approach. In the problems that we deal with, the resulting operator 
is not given but instead the resolvent is renormalized. 

To answer the existence we use the following theorem taken from semi-group theory. Let A 
be a subset of the complex plane. A family J{E), e A of bounded linear operators on the 
Hilbert space "H under consideration, which satisfies the resolvent identity 

J(Ei) - J{E,) = {E, - E2)J{E{)J{E2) (1) 



for £"1 , £"2 ^ A is called a pseudo resolvent on A 

The following corollary (Corollary 9.5 in [12]) gives the condition for which there exists a 
densely defined closed linear operator A such that J{E) is the resolvent family of A: Let A be 
an unbounded subset of C and J{E) be a pseudo resolvent on A. If there is a sequence E'fc e A 
such that \Ek\ — ?■ 00 as A; — )■ 00 and 

lim -EkJ{Ek)x = x, (2) 

for all X G 7/, then J{E) is the resolvent of a unique densely defined closed operator A. As 
we will see, the family satisfies J{Ey = J{E*) so it is a holomorphic family of type (A) in the 
sense of Kato [T3j. Hence, it defines a self-adjoint operator. 

Let us mention the possibility of using the results from [3] in the case of Dirac-delta functions. 
In the approach of [3] we consider an operator A with a dense domain, and consider the singular 
perturbation by an element (p in some dual space, formally: 

A- =A + X{4>,.)4>, (3) 

here the bracket refers to dual pairing. The interesting case is when we have (p G 'H-2{A), 
where 

^en^,{A) if IIy^<^II<^ (4) 

and ||.|| refers to the usual norm in the Hilbert space. Then the theorem in [3] states that the 
operator A~ provides a self-adjoint extension with a new domain. In our case, 

\\————^Sg{a,.)\\= ds s e ' Ksia,a;g) < 00 (5) 

thus we have the same type of singular perturbation-so called form unbounded one. It is 
interesting to see how the conditions on Ricci curvature we found will arise in this approach. 

In our presentation, we do not follow a formal writing style since the paper is rather long 
and has many technical details, hopefully this makes reading more enjoyable. Various operator 
identities that we use can be proved rigorously in the compact case, but they require some 
more work in the noncompact case and it can be done by using the spectral theorem for the 
Laplace-Beltrami operator. We refrain from completing these arguments since they are more 
or less standard in operator theory. 
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2 Point Interactions in Two and Three Dimensional Rie- 
mannian Manifolds 



We adopt the natural units h = 1 in the non-relativistic models discussed in this paper for 
simplicity. In [8], after the renormalization we have found the resolvent kernel corresponding 
to the Hamiltonian for the point interactions (Dirac-delta interactions) in two and three 
dimensional Riemannian manifolds as 

TV 

Rix,y\E) = Ro{x,y\E) + ^ Roix,ai\E)<^r^{E)R^{aj,y\E) , (6) 

where 

{[ dtKt{a,,ai;g) (e-'^" -e'"") if 2 = j 
[~ J Kt{ai, aj; g) e* ii i j . 

Here "^{E) < and Kt{x, y; g) is the heat kernel on the Riemannian manifold, which is defined 
as the fundamental solution to the heat equation 

1 ^2z^/ ^ 9Kt{x,y;g) , . 

—VMx,y;g) = . (8) 

In order to show that the resolvent kernel given in the equation ([6]) corresponds to a unique 
densely defined closed operator H, we need to first prove that it satisfies the resolvent identity, 
i.e, 

R{x,y\Ei)-Rix,y\E2) = iEi-E2) [ z R^x, z\Ei)Riz,y\E2) . (9) 

Jm 

A detailed proof as well as all properties of the heat kernel that we use in this paper, are given 
in our previous work [S] and the relevant literature is also given there. Here we will just give 
the main idea of the proof for the completeness of this paper. If we substitute the equation ([6]) 
into the equation iQ, we obtain 
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Roix,y\Ei) - Ro{x,y 


E2) + Roix,ai\Ei 








N 










= {E, - E2) [ dfz 


Ro{x,z\Ei)Ro{z,y\E2 


Jm 





TV 



+ Ro{x,z\Ei)Ro{z,ai\E2)^^\E2)Ro{aj,y\E2) 

TV 

+ Ro{x,ai\Ei)<^-\E,)Ro{aj,z\Ei)R^{z,y\E2) 



TV TV 

+ ^ ^ Ro{x,ai\Ei)<i>-/{Ei)Ro{a,,z\Ei 

i,j=l r,l=l 



xRq{z, ar\E2)^^i\E2)Roiai, y\E2) 



(10) 



The term Ro{x, y\Ei) — RQ{x, y\E2) equals to the first term in the right hand side of the equation 
above since the free resolvent kernel RQ{x,y\E) must satisfy the resolvent identity (j9]). If we 
add and subtract the terms 
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J2 Ro{x, a,\Ei)<^i/{Ei)Ro{aj, y\E2 



N 



Ro{x,a,\E^)<^;^\E2)R^{a„y\E2) 



(11) 



(12) 



to the remaining terms in the equation above and rearrange, one can complete the proof for the 
resolvent identity by showing that the difference of the principal matrix $jj(i?2) — $jj(-E'i) 
equals to the difference of free resolvent kernels Ro{ai, aj\Ei) — Ro{ai, aj\E2). It is easy to show 
this by using the formula expressing the free resolvent kernel as a Laplace transformation of 
the heat kernel and semigroup property of the heat kernel following a change of variable for 
the time variable in the heat kernel [8]. The equation (j2]) requires the following condition to 
complete the second part of the proof 



\EkR{Ek)f + f\\^0 



(13) 



as /c — >■ oo, where / belongs to the Hilbert space T-L = L^(A^) and the norm is the usual L^(A^) 
norm. Let us choose the sequence E^ = —k\Eo\, where Eq is below the lower bound E^ on the 
ground state energy which has been found in [8]. Then, we must show that 



as — 7- oo. Using the equation (|( 

\\\Ek\RiEk)f 



\\\Ek\R{Ek)f-f\\^0, 

and separating the free part, we get 

/II < \\\E,\R,{E,)f - f\\ 

+ \Ek\\\Ro{Ek)^~\Ek)Ro{Ek)f\\ . 



(14) 



(15) 



It is well known that the first part of the sum converges to zero as k — oo, that is, the 
free resolvent corresponds to a densely defined closed operator (Laplacian). Moreover, the 
Laplacian on geodesically complete Riemannian manifolds is essentially self-adjoint in L^(A^) 
p3|[T5]. Therefore, we are going to investigate only the second term in two and three dimensions 
separately. Two dimensional analysis has been already worked out in [S] and we will just 
review it here and then give the detailed proof for the three dimensional case. Since the 
norm of an operator is smaller than its Hilbert-Schmidt norm: PI I < Tt^^^{A^A) with A = 
Ro{Ek)^-\Ek)RoiEk), we have 

\Ek\\\Ro{Ek)<^-\Ek)Ro{Ek)f\\ 

N 



< \Ek 



/ d^x Ro{ai,x\Ek)Ro{x,ai\Ek) 

dfy Roia„y\Ek)Ro{y,ar\Ek)\^T.\Ek)\ |$,V(^fc) 



1 1/2 



.(16) 
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Let us first consider the diagonal case / = i and r = j for the terms inside the bracket above. 

iV 



\E,\ 



/ dfxRo{ai,x\Ek)Ro{x,ai\Ek) 



d^2/ Ro{a,,y\Ek)Ro{y,a^\Ek)\%^{Ek)\ m\Ek 



M 



< \E. 



1/2 



1/2 



A^^ max ai{Ek) max aji^E^) max 



l<i<Ar 



(17) 



where we have defined Oii{Ek) = fj^d^y Ro{cii,y\Ek)RQ{y,ai\Ek) for simplicity. It is easy to 
see that 



dfx Ro{ai,x\Ek)Ro{x,ai\Ek 



M 



^0 

oo 



dh dt2 A'i,+i,(a„az;^7)e"(*^+*^)l^'=l 
dtti^i(ai,az;^)e-*l^'=l (18) 



by using the fact that the free resolvent kernel is just the Laplace transform of the heat kernel. 
The upper bound of the heat kernel was given in [T6l [T7] and summarized in [8] for compact 
(with bounded Ricci) and Cartan-Hadamard manifolds [9]. We shall use the notation for the 
dimensionless constants coming from the bounds of the heat kernel as C with subscripts for 
simplicity since the exact form of these constants do not play any role here. The upper bound 
of the heat kernel for compact (with bounded Ricci) and Cartan-Hadamard manifolds is given 
in the following form 



Kt{x,y]g) < 



Ci 



V{M) (t/2m)^/2 

C4 / md 

■exp 



(t/2m)«/2 



TfidPix 

exp ( — — ^ — ) for compact manifolds 



(19) 



for Cartan-Hadamard manifolds 



where V{A4.) is the volume of the manifold and d{x, y) is the geodesic distance between the 
point X and y. Then, on-diagonal upper bound of the equation (ITSl) for compact manifolds 
(with bounded Ricci) becomes 



max ai{Ek) < . ... „ 



D , , D 



C,{2m)- \Ek 



4^-2 



where Cg = C2T{2 — D/2). For Cartan-Hadamard manifolds, we get 



max aAEk) < C4 (2m) 2 \Ek\ 

l<i<N 



(20) 



(21) 



We have also 



N 



max |$ri|2 < max ^ l^'.^l^ = max {<^-\Ek)<^-\Ek)h < Pi^'^iEk)) 

l<t,j<N ■' l<i<N — ^ \<i<N 

< \\^-\Ek)\\<\\^-\Ek)\\^ (22) 
where we have used ^\Ek) = ^{E^) for E^ G M and p is the spectral radius. 



In order to find the upper bound for tlie norm of tlie inverse principal matrix, we first 
decompose the principal matrix into two positive matrices 



^ = D-K 



(23) 



where D and K stand for the on-diagonal and the off-diagonal part of the principal matrix, 
respectively. Then, it is easy to see $ = D{1 —D~^K). The principal matrix is invertible if and 
only if (1 — D^^K), and (1 — D~^K) has an inverse if the matrix norm satisfies ||-D~^A'|| < 1. 
Then, we can write the inverse of $ as a geometric series 



$-1 = {I- D-^K)-^D-^ 

= {l + {D~^K) + {D'^Kf + ..)D-\ 



(24) 



and the norm has the following upper bound 



= 11(1 - D~^K)-^D-^\\ < 11(1 - D-^K)-^\\ \\D-^\ 
- 1 - \ \D-^K\\ " " 



(25) 



Since we are not concerned with the sharp bounds on the norm of $ ^ for this problem, we can 
choose \Ek\ sufficiently large such that H/^^-'^-ft'll < 1/2 without loss of generality and get 



\\^-\Ek)\\<2\\D-\Eu) 
Whenever = diag($j"]^, $^2^, . . . , ^^^j^), then 



max |<l>./| . 

l<i<Ar 



(26) 



(27) 



The lower bound of the diagonal principal matrix for compact and Cartan-Hadamard manifolds, 
which was given in [8], leads to the upper bound of the inverse of the diagonal part of the 
principal matrix. Hence, we find 



C7(2m)-Mn~i ifD = 2 

max \^Z^(Ek)\ < X 
i<i<N I - ^ '^^l - 1 C8(2m)-3/2 



\E,\'/'- 



if D = 3 , 



(28) 



for compact manifolds and 



max \<^^\Ek)\ < < 



c,(2™)-M„-(^ 



Cio(2m) 



-3/2 



1/2 



if D = 2 



if D = 3 



(29) 



for Cartan-Hadamard manifolds. Here ^ is a positive constant and defined in [8]. 

If we substitute the results (^U^, and (^U^ into ior D = 2, and take the limit 
k 00, the result goes to zero. Since the norm is always positive, we prove 



\E,\RiEk)f-f\\-^0 



(30) 



as — oo. For the off-diagonal terms, we do not have to make a separate detailed analysis 
since one can easily show that these terms are essentially exponentially suppressed by the factor 
^-y/2m\Ek\d{ai,aj) ^^iq upper bounds of the modified Bessel functions which are given in [Hj. 

Therefore, all off-diagonal terms exponentially vanish when we take the limit A; ^ oo, which is 
enough for our purposes. 

On the other hand, this proof does not work in the three dimensional case as one can easily 
see. In three dimensions, estimating the operator norm by the Hilbert-Schmidt norm is not a 
good way. Instead we will return to the second term in (fT5l) . and show that 



\Ei, 



[ dlx V R^{x,ai\Ek)<l>^\Ek) [ d^z Ro{aj,z\Ek)f 



xRo{x,ar\Ek)<l>^t^{Ek) / d^y Ro{auy\Ek)f{y) 



M 



1/2 



(31) 



goes to zero as — > oo for any / G L^(A^). From our previous argument, we know that the 
inverse of the principal matrix $ satisfies for three dimensional compact and Cartan-Hadamard 
manifolds: 

max \%'{E,)\<^^^^^^^, (32) 

where we define all the constant terms coming from the bounds of the heat kernel as Cu and 
ignore the term in the denominator for large values of \Ek\ for simplicity. Moreover, we can 
combine the two resolvents with the common variable x in the equation (I3T!) . As a result, 
we can express this combination as in the equation (fT8|) and the diagonal upper bounds of it 
for = 3 has been given in the equations (!20l) and (12T]) for compact and Cartan-Hadamard 
manifolds, respectively. Once again, we skip the detailed calculations for the off-diagonal terms 
{i 7^ r) in the above sum since they are exponentially suppressed by the factor Q-y/^MEk\d(ai,ar) ^ 
We always concentrate on the least convergent part in the terms. Once we have achieved our 
goal for these terms, we are done. 

Therefore, it is sufficient to deal with only the diagonal term (r = i) in the equation (13T]) . 
It is easy to show that it is smaller than the following term 



1 1/2 r 



max aAEu] 

l<i<N 



1 1/2 



max |$i/(^fc)| max \^ii\E, 



l<i,l<N 



X 



max 

i<i<JV 



j^dlzRo{aj,z\Ek)\f{z)\] max ( j^d^y Ro{y,ai\E,)\f{y)\ 



1 1/2 



■ (33) 



Using the equations ( l20l) . ( 1211) and ( l32l) in the above equation, we get the upper bound of ([33 
for three dimensional compact manifolds 



N^/^\Ek 



V{M)\Eu\'' 



ISa.P/2 



max 




1/2 



Cn(2m)-3/2 



I^.P/2 



dgZ RQ{aj,z\Ek)\f{z)\ ) max 




dlyRo{y,ai\Ek)\f{y)\ 



1/2 



(34) 
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and for three dimensional Cartan-Hadamard manifolds 



"C4(2m)i" 


1/2 


"Cn(2m)-3/2 






\Ek\'/' 



max 



j^d^gZ Ro{aj,z\Ek)\f{z)\^ max (^j^dlyRo{y,ai\Ek)\f{y)\ 



1 1/2 



(35) 



All these imply that the term 



dlyRo{a„y\Ek)\f{y)\ 



(36) 



M 



must decay at least faster than \Ek\ We now recall that the free resolvent kernel is just 
the Laplace transform of the heat kernel 



POO 

Roia„y\Ek)= / dt e-*l^'=lKi(a,-, (?) 
Jo 



(37) 



so that we can find an upper bound for it by using the equation (fT9|) and evaluating the integrals 
over t 



. . mCi2 

Ro{aj,y\Ek) < — rexp 

d{aj,y) 



Ci3d{aj,y)y/m 
' V{M),/\E^\ 



m(f{aj,y)\Ek\ 



l/2n 



md?{aj,y)\Ek\ 



1/2- 



exp 



-21 



m(f{aj,y)\Ek\ 



l/2n 



(38) 



for three dimensional compact manifolds and 



mCi4 

RQ[aj,y\Ek) < 37— r^^P 



d{aj,y) 



md'^{aj,y)\Ek 
C, 



1/2- 



(39) 



for three dimensional Cartan-Hadamard manifolds. Here we have used the upper bound of the 
modified Bessel function Ki{x) given in [8]. 

For simplicity, let us first consider the generic term which is common for both compact and 
Cartan-Hadamard manifolds and keep the inverse volume term aside for the moment. Then, 
we find for the generic term 



d^y Ro{aj,y\Ek)\f{y)\ < mCu / d^y exp 

M JM 



We now divide the integration region into two pieces 



md\a^,y)\E,\\'/^] \f{y)\ 



C, 



(40) 



dly exp 



-2 



md\a,,y)\E,\\''^^ \f{y)\ 



J d{aj,y) 



M\Bs{aj) 



dgy exp 



md'^{aj,y)\Ek 



1/2 



\f{y)\ 

d{aj,y) 



(41) 
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where Bs{aj) is the geodesic ball of radius 6 centered at aj. It is easily seen that 



M\Bs{aj) 



< -exp 





< - exp 





<-exp 



dly exp 



-2 



m(f{aj,y)\Ek 



1/2. 



\f{y)\ 



I I 



1/2. 



1/2. 



1/2. 



d^y exp 



>!\B,(a,) 



mrf2(aj,?/)|Efc| 



1/2 



\f{y)\ 



dlv exp 



m(P{aj,y)\Ek 



1/2. 



d„?/ exp 



-2 



md^(aj,y)|Efc| 



1/2 



1/2 



11/11, (42) 



where we have used the fact that d{aj, > ^ for all j and y G \-B5(aj) in the second line. We 
then find an upper bound in terms of the norm of the function f{x) by using Cauchy-Schwartz 
inequality in the last line. 

For compact manifolds, it is a simple matter to find the upper bound to the above integral 



M\Bs(aj) 

<iexp 



< - exp 




^ly exp 



md\a^,y)\E,\V'\ \f{y)\ 



C, 



d{aj,y) 



m6'^\Ek\ 
m5'^\Ek\ 

Co 



1/2. 



1/2. 



V{M) sup exp 

y&M \ 

V''\M)\\f\\. 



md^{aj,y)\Ek 
C, 



1/2. 



1/2 



(43) 



due to the fact that the volume of a compact manifold is finite. For non-compact manifolds, 
it is useful to consider the above integral in the Riemann normal coordinates near one of the 
centers Oj. We further assume that the radius of the ball 6 is less than the injectivity radius 
inj(aj). Let us recall that in Gaussian spherical coordinates, the volume integral of a function 
h on a. D dimensional Riemannian manifold Ai can be written as 



[ dfxh{x) 
Jm 



pn 



dr r''-^J{r,e)h{r,e) . 



(44) 



Here 9 = {Oi, . . . Od-i) denotes the direction in the tangent space around a point that we choose, 
and refers to distance to the cut locus of the point in the direction 6. Hence, we get: 



M 



dly exp 



md'^{aj,y)\Ek\ 



C, 



1/2 



pn 



dVt I dr r J{r, 6) exp 



mr^ I Ek I 



1/2 



(45) 



To proceed further we assume that M. has Ricci tensor bounded from below by Ki, i.e. 
Ric(.,.) > Ki g{.,.) everywhere and the sectional curvature is bounded from above by K2 
on Bs{ai). The upper bound on the sectional curvature is automatically satisfied, since there 
are a finite number of Dirac-delta centers and because we take the metric to be C°°{Ai). Had 
we considered a random distribution of Dirac-delta function centers, in which case they could 
have been located at arbitrarily distant points where the sectional curvature could have been 



unbounded, we would have had to constrain the sectional curvature from above. Then Bishop- 
Gunther volume comparison theorems state that the Jacobian factor of the Gaussian spherical 
coordinates in D dimensions satisfies [181 [19], 



sn^-Vr) sn^^Vr) 
^^<J(r,^)<-^, 



where 



sn^(r 



r 



2m(^ ifir>0 

r if = 



(46) 



(47) 



Then the upper bound of the equation fl45l) becomes 

dl] 



°° , sinh2(v/^:^r) 
dr- 



exp 



< 



|Ki|3/2 



vr 



dfi 



,2 
oo 



dr' sinh (r') exp 



dr' exp 



-2r' 



m 









21 7-- l\ 1/2 

mr \Lk \ \ 1 

mr'^ I Ek I 
Csli^il 

- 1 



2 

1/2 



TT 



2|i^i|3/2 



-m|gfc| 
Cal^-il 



1/2 



(4J 



as long as 



1/2 



> 1. Since we are interested in the limit /c — )■ oo it is satisfied for 
sufficiently large values of \Ek\. Therefore, equation ( l42i) is smaller than 

1/2 

'm5'^\Eu 



5\Ki\^/^ 



exp 



1/2 



m 


Ek\ 







1/2 



1 



-1/2 



(49) 



If we choose 6 = {m^/R)^^^^\Ek\~^^^ , where R is an appropriate scale coming from the Ricci 
tensor around a point, where Ricci tensor is non-zero. The prefactor multiplying the exponent 
goes to infinity whereas the exponent decays rapidly. In fact, it decays fast enough to make 
the expression as a whole go to zero as /c — )■ oo. 

Let us go back to the first integral in the equation ( 14T1) and write it in the Gaussian spherical 
coordinates: 



dfi / dr r J(r, 6) exp 



(50) 



Let us now make the observation that there are constants A^, A_, which depend only on 6 and 
Kj^s such that, 



(51) 



for r e [0,5]. For this part of the integral we use the following characterization of essential 
supremum: let us define 



A(e) = ^({r G [0,5] I |r3/2F(r)| > e}). 



(52) 



in 



where /i is the standard Lebesgue measure. Then we have 



Essup \r^/^F{r) \ = inf{e|A(e) = 0} . 

re [0,5] ^ 



(53) 



Let us use now F(r) = J^a dfi|/(r, 9)\, and using Bishop- Gunther bound for the first part as, 

exp 



dr r=^/2 / dl] 1/(7 



_2 I' mr^m ^ 



sn^^ [r) 



,1/2 



(54) 



which is smaller than; 



rS p exp 
Al{K,,Q) / drr3/2 / dfi|/(r,^)| 



_2 mrflSfcl^ 



,1/2 



, . f pS exp 
< A^(/s:i,0)(Essup|r=^/2F(r)|) / dr 

^ re [0,5] ^ Y io 

Al{K\,0) 



_2 ( iinfL^ j 



1/2 



,1/2 



< ( Essup \r^/^F{r)\ 



re [0,(5] 



2(m/C3)V4|E,|i/4 



(55) 



If we take the limit 6 = (my/R) ^^^\Ek\ — > 0, we claim that the essential-suppremum goes 
to zero. To see this, we observe by Markov inequality [20] that 



A(e) < - [ dr |r3/2F( 
^ Jo 

. r „A iV2 



r)| 



< 



dr r 



dr rM / dl] \f{r,9)\ 



1/2 



< - 



1^ 

e 72 



dr 



-sn 



sn^^(^rj 



^^(r) / dfi |/(r, 

§2 



1 1/2 



dl] 



< ^^(47r)VM+(0,i^2) 



< 



1_S_ 
e x/2 



(47r)VM+(0,i^2 



I dr sn|^2(r) / dfi |/(r, 

J§2 



dr r^ / dl] J(r,e)|/(r, 

i§2 



1 1/2 



1/2 



(56) 



For any e > 0, our choice of 6 implies that we can make A(e) = 0. Thus, the infimum goes to 
zero in the limit as /c — )■ oo. As a result we see that the equation ( |3T1) is smaller than 



Cl5 



exp 








(s\Kr\yA 





y/2 


\fm\E,\Y' 1 


-1/2 

11/11 




[\C3\K1\J \ 





(57) 
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and it goes to zero as — )■ oo. The repeated application of the same analysis leads us to 
the same conclusion for the other terms coming from the inverse volume term which has been 
omitted for simplicity. Indeed, all these terms decay with \Ek\ faster than the result that we 
have obtained above. This completes the proof of the existence of the Hamiltonian for point 
interactions in three dimensional Riemannian manifolds. 



3 Relativistic Point Interactions on Two Dimensional 
Riemannian Manifolds 

The resolvent for this system have been found in [9] and it is given by 

R{E) = Ro{E) + R^{E)h'^^"\E)hRo{E) , (58) 



where 



and 



TV 



6t = ^0(-)(a,)x. (59) 



i=l 



1 ^ 

1 POO POO 

Here 0*-~^(x) is the positive frequency part of the bosonic field and Oj stands for the position 
of one of the N Dirac-delta function potential centers and /ij is the bound state energy for the 
single delta center at a^. The operator Xj, called angel operator, was first introduced for this 
purpose by Rajeev in [11] and it obeys orthofermionic algebra. Kt{x, y; g) is the heat kernel on 
the Riemannian manifold, which is defined as the fundamental solution to the heat equation 



^lUx,y-g) = ^^^^^^^. (61) 



dt 

We would like to first check whether R{E) satisfies the resolvent identity ([I]). If we put the 
form of the resolvent in second quantized form. 



R{E) = Ro{E) + Ro{E)b'^<l>-\E)bRo{E) (62) 

into the above resolvent identity dH) and we simplify by purely algebraic operations, to arrive 
at the following identity, 

<^i,{Ei) - $,,(^2) + k{Ro{E,) - Ro{E2))b] = , (63) 

where we stripped off the angels and wrote everything in terms of explicit matrix indices, and 
thus ^{E) = ^ij{E)x\xj ^i^d also hi = (j)^~^^{ai) and similarly for bj. Let us now verify the 
above identity, we note that 

-1 roo poo 

^,,iE^) - $,,(^2) = ^ / ds e-'"l^ \ du [e^^^v^^ - e'''^^\e-^'^" K^{a,, ay, g) . (64) 



19 



Let us work out the other term, acting on no particle Fock space, this is the same calculation 
we have done for the renormalized term. For simplicity, we present the calculation in a formal 
eigenfunction expansion of the Laplace operator (which is rigorously valid for only compact 
manifolds) 



k{Ro{E,)^R^{E2))b] = J2f:i 



1 



1 



Ua- IHq — El Ho — E2 



'-Ho + iox — El Ho + ux — E2 



cr A 

1 



6. 



crA 



$^/IK)/A(«i) 



t^A 



t^A V^o 
1 



WAV Wo 



1 



oJx — El uJx — E2 



1 



-OJx — El OJx — E2 







dss dC 







2v^ 



^0 

°° du 







3/2 



U 



00 /•! /'oo 

ds S / dC ^ / 

Jo 2^71" Jo 



ds^ du e-''/^-'^'''K^{a^, g) 



(65) 



After calculating the C integral, we performed an integration by parts over the variable s. Hence 
we have found the required result f l63l) . 

Similar to the previous problem, let us choose the sequence E\^ = —k\Eo\ = —\Ek\, where 
Eq is sufficiently below the lower bound -E* on the ground state energy which has been found 
in [5] and negative. We now want to show fl5Ul) . Substituting the resolvent equation, written 
in the second quantized language, in this expression we obtain the following: 



lim \\Ek[{Ho - EkY^ + (i/o - Eky^<p'^~\ai)%MEk)<p'^+\aj){Ho - Ek)-']f 



fc— >oo 



The free resolvent, that is the first term in the above equation, 

lim \\EkMEk)f + f\\= lim \\EkiHo - Ej,)-' f + f\ \ = 

fc— >oo fe— >oo 



. (66) 



(67) 



already satisfies the resolvent equation hence we should look only into the second part. To see 
this, note that by the triangle inequality 



hm \\EkRiEk)f + f\\<\im {\\EkiHo - Ek)-' f + f\\ 

- fc— >oo 



fc— >oo 



\Ek[iHo — Ek) 



,){Ho-Ek)-']f\\) 



< lim \\Ek[{Ho-Ek] 



^mHEk)<^>^'■\a,){Ho - Ekr']f\\ . 



(6^ 



We choose a one-particle wave function of the form given below. Even though, this is the most 
general one-particle wave function one can write down, it does not include multi-particle wave 



functions. However, due to the mutually non-interacting nature of the particles involved, the 
total Hamiltonian appearing in the resolvent will be a sum of n identical, individual Hamilto- 
nians in the case of a n-particle state and therefore will decay faster than in the one-particle 
case. 



M 



dJxV^(x)0(")(x)|O) = 5^^(a)4L|O). 



(69) 



A direct computation now reveals that, 

cr 

We verify that the limit 



cr 



(70) 



N N 



Urn ||5^5^(i/o+ |i?;t|)-V^-n«.)*.^-'(^fc)</>^+n«,)(^o + |i?fc|)-' 
i=i j=i 

N N 



< lim |Efc|||VV|$r.i(Efc)| \\{Ho+\Ek\r'ct>'^-\ai)ct>^+\a,){Ho + \Ek\Y^ 
converges to zero. An explicit computation reveals that 



(71) 
(72) 



{a,){Ho + \Ek\r'\^) = J2 



■|o). 



{uj^ + \Ek\)uj^ 

The action of {Hq + \Ek\)~^(p^'\ai) onto this expression leads to 



(73) 



(ifo+ii?.i)-y-H«.)E 



[uj„ + \Ek\)uJ„ 



10) =5: 



We have now a term like \(^^MEk)\ \\F{ai,aj\Ek)\\, which satisfies. 



N 



Y,\%\Ek)\\\F{a,,a,\E,)\\\\< 



Tr|$-i(E, 



1/2 



Trlim) 



1/2 



(74) 



(75) 



This can be used in the above norm, and we get after one more use of the Cauchy-Schwartz 
inequality. 



lim \Eu\\\{Ho + \Eu\)-^cl)^-\a{)%\Ek)(t)^+\aj){Ho + \Ek\)-^ 



fc— >oo 



<N\Ek\ max \%\E, 



l<i,j<N 



N 



EE 



E 



r 



1/2 



(76) 



where we have used 



l/a(ai)^(^)l ^ 



[E 



1/2 



E 



UJ-r 



1/2 



(77) 
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We recall that by choosing k sufficiently large we can make the off-diagonal elements as small 
as we like, while the diagonal elements increase. Therefore, without repeating the arguments 
of the previous section for sufficiently large values of \Ek\, we can show that 



max ^ \<l>iMEk)\ < mMEk)\\ < 2 max \<!>t.\E, 



l<i<Af 



where 



max \^u\Ek)\ < 

l<i<N ' " ^ 



c 



16 



(78) 
(79) 



ln(|Efc|/(m-/i™")) ■ 

and the constant Ciq depends on the class of manifolds under consideration. In the above 
equation m should be superseded on Cartan-Hadamard manifolds by its counterpart mcH, as 
defined in reference [9]. Listed below are the values of this constant for compact Riemannian 
manifolds with positive Ricci curvature, for flat space and for Cartan-Hadamard manifolds. 

27r for fiat and compact manifolds 



C 
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for Cartan-Hadamard manifolds 

c{5) 



Thus, essentially we are faced with the sum/integral: 



We will work this out: First we recall that 

1 r(3) 



u,{u, + \E,\y r(2)r(i) 



c 







[u^ + C\E,\f ■ 



Let us now use the exponential form for the integrand; 



r(3) 



dC 



C 



r(3) 



r(2)r(i)7o ' [oo^ + C\E,\]^ 2r(2)r(i) 

and use subordination for oj^, 



rl poo 

/ dC C / ds s'e-^'-e-^^l^^l 
Jo Jo 



2v^ 



u3/2 



s'^ / 4u—m? u ^— X{u)u 



0) 



(81) 



^2) 



^3) 



(84) 



where A(cr) is the eigenvalue of the Laplacian defined in [9]. If we combine the last exponential 
with |/o.(a)p terms, we get the heat kernel at the same points, Ku{a, a; g) and collecting them, 
we find 



r(3) 



< 



40Fr(2)r(i) 

r(3) 
4v/^r(2)r(i) 

r(3) 
4v/^r(2)r(i) 
r(3) 

4v/^|Efc|2r(2)r(i) 



°° dn 
°° du 



~ du 



rl roc 

dec dss^ e-^'/^"-™'"i^„(a, a; (?)e-^^l^^l 
Jo Jo 

poo rl 

/ ds e-^'/^"-"^'"i^,(a,a;^) / dC Ce 
Jo Jo 



u 



3/2 



du 



u 



3/2 



s'\Ek\^ 

Cl7 , Cis 



\Ek\ 
1 

\Ek\ 



du 
ds s 



C 
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+ 



G 



20 



.A{M) u J 

) 

Cigm C2om^ 



\E,\e-^\^^\ 


poo 

/ dv 




Jo 



,3/2 



1 _e-«|i?/=l - s\Ek\e-'\^'\ 

(ms)^ /Av—v 



,5/2 



5) 



1.^ 



where A{Ai) is the area of the manifold. Here the most divergent contribution comes from the 
last term in the above expression, so we first analyze this term. By inspecting the following 
integral representation of the modified Bessel function K3/2{v) [21], 



1 /ms\ 3/2 



2V 2 / 



dv 



,3/2+1 



-{ms)'^ /iv—v 



(86) 



we obtain the following 



d^ 

7^ 



K3/2{ms) 



(87) 



We now use another integral representation of the modified Bessel function ii'3/2(a;) [21] : 

r(2)23/2s3/2 cos(mt) 



As a result we see that 



< 



< 



< 



\Ek\~ Jo 
C22 



ds s dr s- 



cos{msr) 

(^2g2 _|_ g2\2 



°° ds 

~2 



C23 



°° ds 

~2 



C23 



1 _ (,-s\Ek\ 


— s 


Ek\e 




"1 _ e-^l^^l 


— s 


Ek\e~ 





dr 



(r^ + 1)^ 



°° ds 



\Ek\ Jo ^ 
We now note that the integral 



1 - e " - se'^ 



ds 



1 - e ^ - se~ 



is actually convergent. The first term instead becomes. 



24 



ds 



< 



< 



ds 



Hence, 



\Ek\ 
C25 

Wjo 
C26 
A{M)\Ek\^ ■ 



1 _ e-^l^'^l - slEJe-^l^^l 



1 - e - se~ 



-ms/\Et,\ 



E 



UjJuJa + \Ek 



< 



a 



27 



a 



26 



\Ek\ A{M)\Ek\' 



is shown. As a result we see that 
— \Ek\R{Ek 



lim 

A;— >oo 



< lim \Ek\ 

k—^oo 



28 



29 



Efc|ln(|Efc|) A(Al)|Efc|3 1n(|Efc|;. 
which proves that our formula defines a densely defined closed operator. 



(89) 



(90) 



(91) 



(92) 



^ , (93) 



Ifi 



4 Non-Relativistic Lee Model in Two and Three Dimen- 
sional Riemannian Manifolds 



4.1 The Lower Bound on the Ground State Energy 

After the renormalization of the model in [7], the principal operator was given explicitly in 
three dimensions. We can similarly extend the calculations given in three dimensions to the 
two dimensional case [10], so that we have 

POO 

= Ho~E + fx + X^ dt Kt{a, a; g) [e-'^^->'^ - e-*(^o+— s)] 
Jo 

poo P 

-AM dt d^xd^yK,{x,a-g)K,{y,a-g)c^l{x)e-'^''-^^^-'^^Uy), (94) 
Jo Jm^ 

where D = 2, 3 and yU. is the experimentally measured bound state energy of the system con- 
sisting of a boson and the attractive fermion at the center. In this section we will restrict E to 
the real axis. In order to give the proof that the energy E is bounded from below, we split the 
principal operator as 

^{E) = K{E) - U{E) , (95) 



such that 
and 



K{E) = Ho-E + fi, (96) 

POO 

U{E) = Ui{E) + U2{E) = -AM dt Kt{a, a; g) [e-'^""-"^ - e-*(^o+m-s)j 

Jo 

+ AM dt dfx dfy K,{x, a; g) K,{y, a; g) 0t (x) e-*(^o+2,„-^) <p^{y) . (97) 
Jo Jm^ 

It follows immediately that K{E) > nm — i? + /i, so it is a positive definite operator from our 
assumption E < nm + fi. Due to the positivity of the heat kernel and since the difference of 
the two exponentials is a positive operator, the first integral term Ui {E) is a negative operator. 
We thus remark that 

U{E) < U2{E) . (98) 

This clearly forces 

$(E) > K{E) - U2{E) , (99) 

or rewriting it as 

$(E) > K{EY'^ (l - U2{E)^ K{Ef/^ , (100) 

where U2{E) = K{E)-^''^ U2{E) K{E)~^/^ and K{E), U2{E) are positive operators (so is 
U2{E)). It must be emphasized that the unique square root of the positive self-adjoint operators 
K{E) are well defined for all real values of E below We will now show that by choosing 
E sufficiently small it is always possible to make the operator (^{E) strictly positive, hence it 
becomes invertible, and has no zeros beyond this particular value of E (in the last section, the 
self-adjointness will be further clarified). Therefore, if we impose 

\\U2{E)\\<1, (101) 
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then the principal operator ^{E) becomes strictly positive. For Cartan-Hadamard manifolds, 
we have obtained in [7] 



r(l-f) 



n 2 



(102) 



Then the strict positivity of the principal operator f llOip implies a lower bound for the ground 
state energy 



Egr > nm + /i — I nCsiX'^m^^'^ 



"~D" 



(103) 



where 



C31 — C'so 



7rr (2)r(i-f)^r(2-f) 
r(i)^r(|-f)2 



(104) 



For the compact manifolds with Ricci curvature bounded from below by —K > 0, we have 
similarly obtained 



\U2m\ < n 



A^r(2) 



V{M)fj. 



+ 



4^/i/2^D/4^i/2r(2 - f )r(i - f ) 



^f\/(Ai)V2r(f-f) 



+ 



A'm^/27rr(2 - f )r(i - f )2 



V2 4 



(nm + /i — E)"^ 2 



so the lower bound of the ground state energy was found 



(105) 



Egr > nm + /i — (nX^C-i2) ^ ^ 



where 



a 



32 



r(2) 
r(i/2)^ 



+ 



4^/i/2^z)/4^i/2r(2 - f )r(i - f ) 



y(7W)^f /if v^(M)i/2r(| - f ) 



A'm^/27rr(2 - f )r(i 



D\2 
4 / 



(106) 



(107) 



Therefore, the lower bounds on the ground state energies for different classes of manifolds (11031) 
and fll06p are of almost the same form up to a constant factor, so the form of the lower bound 
has a general character. 



4.2 Existence of the Hamiltonian for the Lee Model in Two and 
Three Dimensional Riemannian Manifolds 

The explicit formula for the resolvent of the Hamiltonian in terms of the inverse of the principal 
operator $~^(_E) is given in O [10] by 



R{E) 



H -E 



a 7 
P 5 



:io8) 
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where 



1 



a 

7 
h 



Hq — E Hq — E 
-^-\E) b ^ 



^ '^-\E) b ^ 



Ho-E 



1 



Ho-E 
<^-\E) 
X(j)g{a) . 



Hq — E 

b^ <I>-\E) 



(109) 



Let us check that the resolvent identity R{Ei) — R{E2) = [Ei — E2)R{Ei)R{E2) is satisfied, 
that is, we must have 



a{E^) - a{E2) 7(^1) - 7(^2) 
/3(Ei)-/3(E2) 5{E,)-5{E2) 

_ ( aiE,)a{E2) + 7(^i)/3(^2) «(^i)7(^2) + 7(^1)^(^2) 

^""1 ""^J I f3(^E,)aiE2) + 6iE^)f3iE2) P{E^h{E,) + S{E,)6{E2) 



(110) 



We first consider the first diagonal element of the above matrix. Using the fact that free 
resolvent satisfies the resolvent identity, we get 



Ro{E^p<!>~\E^) 



$(^1) - $(^2) + b{Ro{E^) - Ro{E2)p 



+ E1-E2 



(111) 



^-\E2)bRo{E2) = 



Let us look at the term in the square bracket more closely. By using the explicit expression of 
the principal operator this term becomes 



PCX) 

/ dt Ktia, a; g) [e-'^H,+^-^,) _ ^-t^n.+m-E,)!^ 
Jo 

POO P 

+ AM dt dfxdfy Ktix,a;g)Ktiy,a;g)(l)l{x) 

Jo J 



^-t{Ho+2m~E2) _ ^-t{Ho+2m-Ei) 



+ AVg(a) 



{Ho - Ei)~' - {Ho - E2)-' 



[a . 



My) 

(112) 



One can shift the operator (pUx) to the left 



Hn-E 



0t(x) 



/ dfx' 4>l{x') r dt e-*(^o+™-^) K,{x, x'- g) , (113) 

J M Jo 



and shift the operator (pg{x) to the right 

/ dfx' rdte-*(^«+"^-^)i^,(x,x';<7)0,(x') 
io — -t^ Jm Jo 



^9[^) rr 

no — 



(114) 



IP 



which we have also used in [7] for the renormahzation. The last term in the equation ( I112p can 
be normal ordered as 



/ dt Kt{a, a- g) [e-^i^o+^'^-^i) _ ^-KHo+m-E,)^ 



JM 



A'/ dt I dfxdfy Kt{x,a-g)Kt{y,a-g)(Pl{x) 



Then we prove that 

- $(^2) + 6(^0(^1) - ME2)P + Ei-E2 = 0. 
The other term in the matrix equality flllOp 



7(^1) - 7(^2) = {El - E2] 



a(Ei)7(E2)- 7(^1)5(^2) 



;il5) 



;ii6) 



;ii7) 



can be written as 



- Ro{E,p^-\E,] 



- $(^2) + 6(^0(^1) - Ro{E2)p + E1-E2 



$-1(^2) = 0, (118) 



due to (11161) . Similarly, the other terms can be put into the following forms 



'^-\E,] 



<^{Ei) - $(^2) + b{Ro{Ei) - Ro{E2)p + E1-E2 
$(Ei) - $(^2) + h{Ro{Ei) - RoiE2)P + E1-E2 



<l>~\E2)bRo{E2) = 



<^-\E2) = 0, (119) 



and they are all satisfied thanks to the equality flll6p . Hence, we prove that the resolvent 
identity is satisfied. 

Recall that the resolvent for the Lee model is defined in the following Fock space J-b^""^^-* {H)® 
X+ © Tb^^\'H) ® X-y for any given n G N, and x± is the spin states. In matrix form, we have 
RiE) : J^t3^''+^\n) © J'b^''\H) J^b^'^^^KH) © J'B^''\n). Then we must show that 



\E^R{E,)\f) + |/)|| = \\\E,\R{E^)\f) - \f)\\ ^ 



(120) 



as /c —7- 00. Here |/) G J^b'^^^^\T^) © ^b^^\T^) and the norm is taken with respect to 
Jb^^+^H^) ©-^6^"^"^)- Let us decompose the vector |/) as 



where 



|j(n+l)^ 
|/{n)) 



d^Xi . . .df*X„ f{Xi,X2, . . .,Xn)\xi,X2, ■ ■ ■ , Xn) ■ 



:i2ii 



(122) 
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So we have 



\Ek\a{Ek) \Ek\^{Ek) 
\Ek\(3{Ek) \Ek\S{Ek) 



|j(n+l)^ 



|j{n+l)^ 



n 1/2 

+ + \E,\5{E,)\f(-^) - 



(123) 



1/2 



since ||y4 + B||<||A|| + ||i?||. We shall investigate each norm separately. Let us first consider 
the term 

\\X\E,\<^-\E,)Ma) \ |/(-+^))|| < X\E,m~\E,)\\\\<j>,{a) \ |/("+^))|| . (124) 



Using the formula f lll4p for = — l-Efcl and x = a, we get 
1 



Il0.(a) 



|/^"+^))||<|| / dfx / dte-*l^^li^,(x,a;^7)0,(x)|/("+^))|| 



< 



d-x 

M \ JO 



dt e-*l^'=lKi(x,a; g) 



1/2 



(125) 



Let us first consider the compact manifolds. Then, one can take the integral over the variable 
t by the help of the upper bound of the heat kernel fll9p for compact manifolds and obtain 



\\^g{a 



X 



M 



Hq + \Ek 



^) (v^^^ (2rf(x,a)v/HE^3 

Ko_^ {2d{x,a).Jm\Ek\/C^ 



\Ek\ J \ 



(126) 



(P{x, a) 



1/2 



We now choose Riemann normal coordinates around the point a, assuming that 5 < inj(a). 
Then, we split the integration region into the two parts as = /g^(^) + Xv(\B^(a)- Expressing 



91 



the first integral in the Gaussian spherical coordinates, we get 



dfi dr r 



C: 



33 



C: 



34 



L>/4 



(127) 



+ 



dfx 



m(P{x, a) 



\Ei, 



-^^Kr(2d{x,a)^/m\E,\/a 



a 



34 



(i^(x, a) 



Kd^^ (2d{x,a)^m\Ek\/Ci 



1/2 



where we have used the equations ( 146|) and (15T]) . Let us now consider the first integral. It is 
smaller than the following expression 



Jo 

I Z7 I \ -0/4 



33 



C: 



34 



V{M) 

2 



:i28) 



One can evaluate the integrals [22] 

dr r^+^fs:i2(ar) 



dr r -R'd^i (a?" 



Ar(i + f)r(2 + f)r(f: 

7r(L)-2) csc(7rL)/2) 
4^2 



(129) 



Jo ^ 



2^r(f) 



' 2 +^ 



{D + 2)a-'- 

where a G M"*" and D = 2,3. Then the upper bound of the first integral in fll27p becomes 



mA^-^(i^i,o)/ a 



35 



\Ek 



[m\Ek\y 



(g+2) 



+ C^eimlEk 



V — -1 



37 



y(7W) 



(130) 



For Cartan-Hadamard manifolds, we do not repeat the analysis above because the upper bound 
of the heat kernel for Cartan-Hadamard manifolds given in the equation ( !T9|) corresponds to 
removing the volume term from the one for the compact manifolds. As a result, we get the 
upper bound of the first term in the equation (11271) for Cartan-Hadamard manifolds 



\Ek\ 



m\E,\)^-^ . 



(131) 



9.9. 



Let us now consider the second term in the equation (11271) for compact and Cartan-Hadamard 
manifolds. Due to the upper bounds of the Bessel functions used in [5], we find for compact 
manifolds 



M\Bs{a) 



dfx 



mdP{x^ a) 



v{My 



34 



m\Ek\ 
(P{x, a) 



D/A 



< 



X 



M\Bs{a) 



dfx 



m(P{x^ a) 



\Ek 



C: 



33 



V{M) 



exp 



Kd^^ \ 2d{x,a)^m\Ek\/Ci 
-d{x,a)^/m\Ek\/Ci 



(132) 



2d{x,a)^m\Eu\/C^ 2 



34 



Tn\Ei,\ 



2 1 2d(x,a) 
exp ( -^4^ 



^m\Ek\/C^y 



(P{x, a)) (2d{x, a)v/m|E,|/C3)(4-^)/2 



Since d{x, a) > 6 for all x G \ Bs{a), the upper bound of the above equation is 



exp( -6^/m\Ek\/Cs) [ dfx 

\ J JMKBsia) 



md'^{x, a] 
\Ek\ 



C: 



33 



X 



26y/m\Ek\/Cs 2 



34 



m Eh \ 



52 
X exp 



V{M) 
2 



exp 



d{x, a) 



'm\Ek\/C^ 



{26^m\Ek\/C,)^^-'^y^ 



d{x, a) 2d{x, a) 
2 (4-L)) 



VMEk\/Cs 



(133) 



For compact manifolds, we have a simplification. This upper bound above is smaller than 



exp - 6y/m\Ek\/C3^ dfx 



■md'^{x, a) 



\E, 



a 



33 



V{M) \26^m\Ek\/Cs 2 



+ Q 



34 



m\Ei,\\^^'^ 



J {26y/m\E,\/CsY^-Dy' 



(134) 



Due to the fact the geodesic distance between any two points on the manifold and the volume 
of the manifold is finite, that is, d{x,a) < dmax{a) = maxxd{x,a), the upper bound to the 
above integral can easily be found as 



\Ei, 



exp{-6^/m\Ek\/C3)ViM) 



33 



1 



^ 1 



ViM) \26^m\Ek\/C3 2 



C: 



34 



mlEi-W^^^ 



-\ 2 



5' J {2s^m\E,\/c,y^-^y^ 



(135) 



For Cartan-Hadamard manifolds, we similarly find 



exp 



M\Bs{a) 



if. 



- l¥g|)0"l-B*|/C5 



d'^{x, a) 



< 



Jm 



d X exp 



a) 



4-D 



1 ^/m\E,\/C, 



(136) 



where we have used d{x,a) > 6 for all x G \ Bs{a). Let us write the above integral in 
Gaussian spherical coordinates as we did in Section 2, 



Pn 



an / dr V(r, ^) exp 
'o 



4-D 



m\Ek\/C5 



(137) 



To proceed further we assume that Ai has Ricci tensor bounded from below by Ki. As a result 
of this and using the equations (H6l) and the upper bound to the equation (I136p becomes 



< 



C4om^-i|Efcp-3exp (^-5^m\Ek\/C^ 

52(_Ki)(^^-i)/2 X 



X exp 



dr sinh^ ^ 
4 



-i^ir) 



4-L) 



1 x/m|E,|/C5 



(138) 



Since sinh ^(x) < e^^ j'^p the integral and get 



;-i^i)(«"i)/252 



- 1) V^li^fcl/Cs - (/^ - l)v^ 



(139) 



as long as 



- 1) ^m\E^\IC, -{D- l)V^i 



> 0. But this is always satisfied for 



sufficiently large values of \Ek\. 

Therefore, if we combine the results fll30p and f ll35p we obtain for compact manifolds that 



Ho + \Ek 



■|/("+^))||<v/;^|||/"+^)) 



mAX-\KuO) C'33 , IP |._(iZ±2) 



\Ek\ 



+ C34(m|Efc|)t-i + -^(m\Ek\)-' 



V{M) 



\E,\ 



X V{M) 



33 



V{M)\2Sy/m\Ek\/C3 2 



34 



I Z7 I \ -0/4 



exp ( S^ym\Ek\/C3 
2 



(2V"^l^fe|/C3)(^-^)/2 



1/2 



(140) 
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and the results fll3ip and fll39p for Cartan-Hadamard manifolds give 



\\(pgia) 



Ho + \Ek 



+ 



-|/("+i))||<v/^|||/("+^)) 



38 



11/2 



-K,)(^~^)/^^ (i^ - 1) V^l^fcl/Cs - P - l)v^ 



(141) 



We are now going to find an upper bound of the inverse norm of the principal operator. In 
order to do this, let us recall that we split the principal operator when we try find the lower 
bound of the ground state energy. We now split the principal operator in the following way: 
$ = {K — Ui) — U2, where Ui and U2 are defined exactly as before. Then, we have 



$-1 = {K - Ui 



-1/2 



1-{K- Ui)-^'^ U2 {K - U: 



-1/2 



{K - U, 



-1/2 



(142) 



Let us substitute the identity operator K^/^K between the operators {K — Ui) ^^"^ and V^- 
Hence, 



$-1 = [K - Ui)-^^^ [1 - X]-^ {K - Ui) 



-1/2 



(143) 



where we have defined X = (K — Ui) ^^"^ K^/'^U2K^^'^ {K — Ui) ^^"^ for simplicity. Here the 
following operator can be written as an infinite geometric sum 



(144) 



as long as ||X|| < 1. This leads to 

Il[l-^]"M|<[1-Il^ll]"' • 
Since —Ui is a positive operator, {K — f/i)~^/^ < K~^/'^. Then, we have 

ll^ll < \\U2\\ ■ 

If we make | \ sufficiently large then \ \U2\ \ < 1/2 and 



[l-ll^l 



1-1 



< 2 



As a result of this, we get 



\^-\Ek)\\<2\\K~''\Ek)\\^ < 



\Ek. 



where we have used 



K-"^{Ek) = {Ho + ^i+\E,\)-^'' < 



1^.11/2 



(145) 



(146) 



(147) 



(148) 



(149) 



Then, substituting the equations fll30p and (11350 for compact or substituting the equations 
(1131 p and fll39p for Cartan-Hadamard manifolds into the equation f ll27p and using the above 



9.^ 



upper bound for the inverse principal operator, and taking the hmit as /c — )■ oo, we eventually 
obtain 



\E,\ \\^-\E,)\\ ||0,(a) 



1 



Hq + \Ek 

Let us consider the other terms in the equation fll23p now: 

\Ek\ 



|/(n+l))||^0 



(150) 



|||E,|a(E,)|/("+i))-|/("+i))||< 
1 



Ho + \Ek\ 



Ho + \E,\ 
\\^^\E,)\\ ||0,(a) 



1 



(151) 



Ho + \Ek 



|;(nH-i))|| 



The upper bound for the norm Il77;^qq^</'g('^)ll can be similarly found, and comes out to be the 
same as the one for ||0g(a)-j^^jqq^||(the norms are different in general). As a result of this, we 
obtain 

|||£;fc|a(Efc)|/('^+i))-|/("+^))||^0, (152) 
as oo. Similarly, the following term 



\\\E,\^iE,)\f-^)\\<X\E,\ 



r0l(«)l/^"^)ll \\'^-\E,)\\ 



and 



'Ho+lEkl'' 

|||i^.|/9(i^.)|/^"+^))|| < A|i^.| ||0,(a) ^ \^ , |/(-+^))|| \\^-\E, 



Ho + \Ek 



(153) 



(154) 



both vanishes as — )■ oo. Moreover, we have 



1/^-^)11 



II (|£;.|<5(^.) - 1) |/W)|| = II [|i5;,| <^-HE,) - 1] |/("))| 

\E,\ K-^/\E,) [l + (1 - U{E,)r^ - l] K-^'\Ek) - 1 

< II {\E,\K-\E,) - 1) |/W)|| + \\\E,\K-y\E,) [(1 - U{E,)Y^ - l] K-^/\E,)\f-))\\ 
= II {\Eu\K-\E,) - 1) |/("))|| + |||E,|K-V2(i^,)(i _ U[E,)r^U{E,)K-y\E,)\f^-))\\ ,(155) 



where we have used the fact that the factor (1 — U{Ek)) ^ can be considered as an infinite 
geometric sum. The first term goes to zero as -t- oo since 



\Ek 



Ho + \Ek\+fx 



< 



< 



\Ek 


+ /i 


Ho + \Ek 


+ /i 



1 J 1/^-^)11 
- 1] 1/^"^ 



\Ek 


+ /i 


Ho + \Ek 





1 1/^-^)11 + 



l-^fcl + /i - /i _ 

Ho+\Ek\+fi 
Ho + \Ek\+fi 

\Ek 



1/^-^)11 



(156) 



where the term containing Ho vanishes as /c — )■ oo because it is the free resolvent and the second 
part clearly goes to zero. This shows that the first term in the equation 01551) vanishes in the 
limit. As for the second term, it is smaller than 

\E,\\\K-'/'{E,)\\ [{l-\\U{Ekm-'] [||f/i(i^^.)|| + ||f/2(i^;.)||] ||i^-^/^(i^;.)|| |||/("))||. (157) 

9fi 



Since m > /i for bound states, one can easily see that 

^0 

X {Ho + m+\E,\)-'/''\\ 
^ ^'C43A2||(iJo + m+|E,|)-Mn(tedM)|| forD = 2 
CiA^\\iHo + m+\Ek\y'/^\\ forD = 3 ' 



(158) 



Here we use the fact that the operator in the parenthesis, which we call A{s) is positive, and 
for a positive family, if two integrable functions satisfy < /(s) < g{s), then J ds f{s)A{s) < 
J ds g{s)A{s). Moreover, for positive operators, order relation implies the same ordering for 
the norms of the operators. For simplicity, we have also disregarded the more convergent in 
\Ek\ coming from the volume terms of upper bound of the heat kernel for compact manifolds. 
We now note that 



In 



Ho + \Ek\+m\ .. ,, dt 



Hq + \Ek\ + m \ m — fi J Jq [Hq + \Ek\ + m]t + m — fi 



< 



< 



< 



dt 



[Hq + \Ek\ + m]t'^ + m - jj, ,^ggs 

1 r ds ^ ' 



^/W^[Ho + \Ek\+mYl'^'' Jo s2 + l 
C45 



^/m^^\Ek\^/'^ ' 



Hence, we get 



cUeX^ for D = 2 

\Ui{E,)\\ < {^^^ „ ^ . (160) 

for = 3 

Using the results f ll02p and f llOSp for E = —\Ek\ with the above analysis, we finally obtain 

\E,\\\K-\E,)\\ \\UiE,)\\ \\il-UiE,)r'\\ |||/("))||-^0, (161) 

as /c — 00 and this completes the proof that our renormalized formula corresponds to the 
resolvent of a densely defined closed operator. 

We will further show that $(-E) is a holomorphic self-adjoint family of type (A) in the sense 
of Kato [13]. This will in turn justify the claim that the resolvent corresponds to a self-adjoint 
operator. To prove this we will use the theorem given by R. Wiist |23]- First, we define a 
holomorphic family of type (A) as follows: Let G* C C be domain and L{z) be a family of closed 
linear operators, acting on a Hilbert space "H, {L{z)\z G G}. If 

1) the domain V{L{z)) = V is independent of z G G 

2) for any / G and g then {g\L{z)\f) is holomorphic in G , (162) 

then this is a holomorphic family of type (A). 

An operator which is a holomorphic family of type (A) is a self-adjoint holomorphic family 
of type (A) if 

1) G is a symmetric domain of the complex plane relative to the real axis. 
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2) V is dense in "H 

3) V{L{z)^) = V{L{z*)) for a\\ z e G . (163) 

Theorem(Wiist): Let G be a symmetric domain of complex plane relative to the real axis, 
and L{z) is a holomorphic family of type (A) defined on G. Assume 

1) V is dense in T-L 

2) V{L{zy) D V{L{z*)) . (164) 

Let 

M = {z e G\L{zY = L{z*)} . (165) 

If M is not the empty set then it is the whole domain G. This implies that the family is a 
self-adjoint holomorphic family of type (A). 

Let us consider our case. We choose the domain G as 

G = {Ee C\^{E) < fi} , (166) 

which is symmetric with respect to the real axis. The principal operator $(£') given explicitly 
in ([MD formally satisfies the relation $(E)t = $(E*) so this implies V^^EY) D V{^{E*)). 
Let us assume that the family is holomorphic for now. Note that a densely defined holomorphic 
family of operators satisfying the formal relation ^{E)^ = ^{E*) is closable. Proof: Let us 
consider the common domain V, and choose \gi) G P ^ as / — )■ oo. We further assume that 
^{E)\gi) converges to some \g{E)). Then we have 

{fmE)gi) = {<!>HE)f\9i) = {HE*)f\gi) ^ , (167) 

for any |/) G I' as / — )■ oo. This implies that 

^E)\gi) ^ \g{E)) = . (168) 

This is the requirement of closure. Of course, we must establish this closure uniformly, that is, 
we need to show that for every sequence \gi) converging to some element, if $(£'o)|5'i) converges 
for one Eq inside the region G, then it converges for all E E G. Hence, we can define a unique 
closure over G, the closures having a common domain V~. Once we determine a common 
domain for the family $(-£'), we will prove that there is indeed a closure over a common 
domain. 

For the family ^{E), we choose V = V{Hq). It is well known that if is a geodesically 
complete manifold then Laplacian defined on is a closed, densely defined self-adjoint operator 
[m [15]. Then, the operator Hq = Jj^ d^x (pl{x){—^V'^)(f)g{x) defined over the direct prod- 
ucts of n copies of the Hilbert space L^(A^) is also a densely defined (essentially) self-adjoint 
operator. Moreover, the finite direct sum of such operators will preserve this property. 

The first term Hq — E + fi is obviously defined over this domain V{Hq) and it is a closed 
operator. Note that the other term can be defined by the spectral theorem, 

oo 

ds Ks{a, a; g) [e-^^™"/^) - e-'*(^o+'"-^)] (169) 



and it is a positive operator when E is real and 9?(£') < yU. Its domain of definition includes 
V{Ho) when ^{E) < fi. To see this we will use a different integral representation. 



POO 

/ ds Ks{a, a; g)[e-'^"'-^'^ - e-^(-f^"-^+'")] 
Jo 
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= ds s Ks{a,a;g) du e-^"(^«+'^-^)-^('""^)(ifo - ^ + /i) • (170) 
Jo Jo 

When the operator acts on an element l/*-"^) in the domain of Hq, the norm of the resulting 
vector is smaller than, 



1 

ds s Ksia,a;g) I dn He-'^^^o+^-^^-'^^-^^IHKi/o - ^ + . (171) 

Jo 

Now we can estimate the following factor using the bounds on the heat kernels given in fll9p . 

"1 



ds sKs{a,a; g) / dM||e 



-su(Ho+ii—E)—s{m—ii) I 



JO 

pi PCX) 

< du ds s 
Jo Jo 



V{M) (s/2m)^/2 



^—sunm^—su(fi—?R{E))—s{m—fi) (^172) 



Thus we show that 



d. Ks{a, a- g) [e-(^o-i^+m) _ ^-si^m-,)^ 1 1 < p^^ _ gfj(^)) (| |ifo|/^"^) 1 1 + |/x - E 



(173) 



where 



Ci / 1 1 

C2(2m)^/2r(2 - f ) 



(nm + /i-3?(E)(f -1 



((n + l)m-3f?(^))T-i-(m-/i)^"i . (174) 



As a result the domain of this operator family includes V{Hq). In fact, by the spectral theorem 
the operators so defined are closed, when we restrict them to a smaller domain, i. e. to V{Ho) 
they remain closed. So the sum of the two pieces, Hq + fi — E and the term above, defined over 
V{Ho) is closed, since they were already closed operators defined over a common domain. 

The last part requires more work, for this we will first show that U (E) is relatively bounded 
with respect to Hq hence its domain includes V{Hq). Moreover, if we have a holomorphic family 
of operators defined over a dense domain, then they are preclosed, that is we can define the 
closure of this family, as we have shown. It is easy to see that 

\\U{E)H,'Ho\f^-^)\\ < \mE)H,'\\ \\Ho\f^-^)\\ , (175) 

where the first norm can be estimated by exactly the same method developed in [7j. So we are 
giving the result in order not to repeat the similar calculations, for n > 1, 

poo 1*1 

\\U{E)H,'\\ < Csn ds dus Kli'{a,a;g)K'J\a,a;g)e-^^"^e'^''^''K (176) 

Jo Jo 

After using the upper bound of the heat kernel given in f|T9l) and defining new variables p = 
Ci/V{M) and q{s) = C2I [s /2m)^'^ we get 



noo n 1 

\\U{E)H^^\ \ < C^sn / ds dus e-'""*e 
Jo Jo 



snm suSti^E) 



p +pq{s) 



yZ/2 ^3/2 



9Q 



1 1/2 




3/4 



(177) 



Taking the s and m integral, we obtain 

nC49 



< 



n 



(2m) ^/^C, 



50 



n(2m)^/^C5i 



V{M){nm-^if {nm-fiy~^/^ v^VXM)(nm - ^) 2-^/4 



:i78) 



since < fi. Thus we choose the domain of U{E) as V{Hq), and now the family is closable 
over this domain. However, as a result of the closure, the domains for different values of E may 
become different. In fact, this does not happen, as we will see. 

Now we show that we can perform the closure uniformly, as a result of the following: for 
any Ei, E2 E G ^{Ei) — $(-^2) becomes a bounded operator. A short computation shows that. 



||$(£;i) -$(^2)11 < 1^1-^2! 



[n 



1)A^ / ds sKs{a, a; g)e' 



-snm^—s{m—fi) 



(179) 



If \gi) G P is convergent to a vector |/), and assume that ^{Ei)\gi) converges to \g{Ei)) for 
one El, then we set $(£'i)|/) = \g{Ei)) to define the closure at point Ei. Then, for any E2, 
we have 



mE2)\gi)- 

= ||[$(E2) 

\\mE2] 



< 



<l>{E2)\f)\\ = mE2)\gi)^ 
- ^Ei)]\gi > +^E,)\gi) 

-<^>(^i)]|||||^?^)- 



- \g{Ei)) - [$(^2) 

- \g{Ei)) - [$(^2) 
^Ei)\gi) - \g{E,) 



<^>(^i)]l/)ll 
^ , 



(180) 



and this shows that whenever l^f;) converges to |/) and ^{Ei)\gi) converges to \g{Ei)), we have 
$(£"2)1(7;) becomes convergent and the resulting vector is exactly equal to $(£2)]/) as it should 
be for the requirements of the closure. Hence the sum of all these three parts will make a 
holomorphic family $(£) with a dense common domain V{Hq). Moreover, the sum is closable 
over a dense common domain which we call V{Ho)~ . 

We would now make holomorphicity more precise, up to now we have not actually made 
use of it. To prove that the family is holomorphic we will refer to the following theorem, which 
is stated in a slightly simplified form according to our needs and the proof of which can be 
found in [21]: Assume X is a measure space with a cr-finite measure i' defined on it, let J be a 
measurable subset of X. Let G be a open domain of the complex plane. Consider a function 
7 : / X G 1-^ C such that 



1) jix,.)eL\x,\i,\) 

2) li-jz) is holomorphic in G 

3) J |di^||7(x, z)\ is bounded on all compact subsets of G 



Then the function 



T{z) = J dh''j{x,z) is holomorphic in G 



:i8i) 



:i82) 



To use the above theorem, let us write our family in the following form = V{E){Hq — 

E + n) = [l + Vi{E) + V2iE)]{Ho -E + fi), where 

POO 

Jo 

poo 

Jo 

By using our previous estimate in (166) and (167) we see that, the first term Vi{E) is indeed 
uniformly bounded in ^E < fi. Hence, the integrals {f^"'^\Vi{E)\g^'^^) for all functions f,g in 
the Fock space are absolutely convergent. Moreover, any such matrix element satisfies all the 
other conditions on holomorphicity and integrability. The second part, again using ideas very 
similar to the previous estimates, can be written as 

/oo rl 
dss du 0j(i^.„(., a; ^))e-^(^°-^+'"+(i-")^+"^"V.(i^s(., a; g)) . (184) 





This integrand as a function of E is holomorphic in E for ^E < fi and it is absolutely integrable 
for any ^E < fi on [0, oo) x [0, 1]. Similarly, it can be shown that the following bound holds 

\{f^-^mE)\g(-^)\ < n- , (185) 

[nm — ^JiE)^/'^ ^ 

which clearly shows that for ^E < /i is uniformly bounded everywhere, hence on compact 
subsets as well. Hence, applying the theorem stated above we see that the resulting function 
V2{E) is holomorphic for ^E < fi. There is one subtle point about the closure operation, 
but this is also solved by the following observation. Let us consider the limit of ^{E)\gi) as 
/ — 7- CX3 as a function of E for any convergent \gi) sequence in the closure operation. If this 
sequence is uniformly convergent on compact subsets the limit is a holomorphic function (by 
an application of Morera's theorem). Note that this family ^{E)\gi) is norm bounded by a 
constant multiple of a simple function given in the equation (171), \Ei — £'2 1- This function 
itself is uniformly bounded on compact sets centered around any given point Ei, hence the 
sequence of functions ^{E)\gi) is a uniformly convergent sequence. This shows that the closure 
remains a holomorphic function of E for ^E < as required. Thus we complete the proof that 
the closure of the family $(£') over D{'Hq)~ is holomorphic for '^{E) < fi. 

Now we are ready to apply the theorem of Wiist. If we choose E' G M and sufficiently small 
E < E^:, then U{E) has relative bound with respect to K{E) which is less than 1. Hence, by 
Kato-Rellich theorem |T3] of perturbations of self-adjoint operators, U (E) will be self- adjoint 
for E < E^. By the theorem of Wiist, the family is self-adjoint everywhere as desired. This 
result is important to establish that the spectrum only lies along the real axis, and justifies our 
search for the lower bound of energy and shows that the resulting operator is self-adjoint as it 
should be. 



5 Conclusion 

In this paper, we have proven that for the three models that we have constructed, namely non- 
relativistic point interactions in two and three dimensional Riemannian manifolds, relativistic 
point interactions in two dimensional Riemannian manifolds and non-relativistic Lee model in 
two and three dimensional Riemannian manifolds, the Hamiltonian after renormalization is a 
densely defined self-adjoint operator. 
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